Critical temperature and Ginzburg-Landau equation for a trapped Fermi gas. 
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We discuss a superfluid phase transition in a trapped neutral-atom Fermi gas. We con- 
sider the case where the critical temperature greatly exceeds the spacing between the trap 
levels and derive the corresponding Ginzburg-Landau equation. The latter turns out to be 
analogous to the equation for the condensate wave function in a trapped Bose gas. The 
analysis of its solution provides us with the value of the critical temperature T c and with 
the spatial and temperature dependence of the order parameter in the vicinity of the phase 
transition point. 

The recent progress in the studies of ultra-cold trapped Bose gases and the discovery of Bosc-Einstein condensation 
||-|| stimulate an interest to macroscopic quantum phenomena in trapped Fermi gases. The most prominent phenom- 
ena should be connected with a superfluid phase transition associated with the appearance of the order parameter 
- macroscopic wave function of strongly correlated two-particle states on the Fermi surface (Cooper pairs). The 
possibilities of finding this phase transition in trapped Fermi gases have been discussed in Refs. 0-6 1. 

A remarkable feature of neutral-atom Fermi gases is that the Cooper pairing and, hence, the superfluid phase 
transition can occur for both attractive and repulsive interparticle interaction. For an attractive interaction (negative 
scattering length a) the pairing occurs in the s-wave channel, as described by the standard BCS approach, and 
one has "singlet" Cooper pairs. In a neutral-atom Fermi gas, according to the Pauli principle, such a pair can 
only be formed by two atoms which arc in different hyperfine states. Therefore, the critical temperature T c of 
the transition is very sensitive to the difference in concentrations of the two hyperfine components, and under the 
i condition An/n = (ni — n-z)/{n\ + ri2) ^ T c jep <C 1, where Ef ~ h 2 n 2 ^ 3 /m is the Fermi energy, there will be a 

complete suppression of the spin-singlet pairing. For 6 Li with a w — 1140A R], one has T c w 30nK for the atom 
■ density n — 4 • 10 12 cm~ 3 [Q, and the existence of the s-wave pairing requires An/n < 10~ 2 . 

For positive scattering length (repulsive interaction) the s-wave pairing is impossible, and one has to consider the 
mechanism of p-wave "triplet" pairing which originates from the effective interaction caused by polarization effects 
^ Actually, this pairing mechanism is insensitive to the sign of the scattering length. It works equally well for 

a < in the situation where the direct s-wave pairing is suppressed. In these cases the Pauli principle allows to have 
t-H ' Cooper pairs formed by two particles which are in one and the same hyperfine state, whereas the particles in other 
hyperfine states participate only in the formation of the effective pairing interaction. Therefore, the p-wave "triplet" 
pairing does not require any severe restriction on An. As found ||, the corresponding critical temperature T c depends 
non-monotonically on n±, ni and becomes zero only in the case where all particles are in the same hyperfine state. 
Since the effective interaction based on polarization effects is weaker than the direct interparticle interaction, the 
p-wave "triplet" pairing results in a lower value of T c as compared to the critical temperature for the s-wave pairing, 
i For example, for the p-wave pairing in a gas of 6 Li the critical temperature T c m 30nK corresponds to densities 
£H 1 n rj 10 13 cm~ 3 . 

In this paper we study the influence of a (harmonic) trapping potential on superfluid pairing. We derive the 
O ■ corresponding Ginzburg-Landau (GL) equation for the order parameter, assuming the critical temperature T c much 
higher than the level spacing Q in the trap. The analysis of this equation provides us with the value of T c for the 
trapped gas and gives the coordinate and temperature dependence of the order parameter in the vicinity of the phase 
transition. As found, the critical temperature is slightly lower than that for a spatially homogeneous Fermi gas with 
density no (maximum density of the trapped gas) , and the behavior of the order parameter resembles the behavior of 
a trapped Bose condensate. 

We consider a two-component neutral gas of fermionic atoms, with a short-range interatomic interaction, trapped 
in a spherically symmetric harmonic potential. The two (hyperfine) components arc labeled by indices a = ± and 
are assumed to have equal concentrations. The Hamiltonian of the system has the form (7i = 1) 



H = /E^W (-^ 2 -/i+^^^a(r) + |^^^(r)^(r)^(r)^(r 



(1) 
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where \x is the chemical potential (Fermi energy), Q the trap frequency, g — Ana/m the interaction strength, and m 
the atom mass. 

The interaction effects can be expressed in terms of a small gaseous parameter A = 2\<i\pf/'k < 1, where pf — 
mvp — (37r 2 no) 1 / 3 is the Fermi momentum. In the spatially homogeneous case the system of fermions described by 

the Hamiltonian (|l|) with il = 0, undergoes a superffuid phase transition. The transition temperature and the 
type of pairing depend on the sign of a. For negative a (attractive interaction) there will be the s-wave "singlet" 
pairing, whereas for positive a the p-wave " triplet" pairing should take place (see Ref. || for details) . In both cases 

the critical temperature — Cep exp{— 1/r}, where ep = p F /2m is the Fermi energy, T the pairing interaction, 
and C a numerical coefficient of order unity. The value of C and the expression for T depend on the type of pairing. 
For the "singlet" pairing one has T = A, and for the "triplet" pairing T w A 2 /13. 

In the Thomas-Fermi approach {ep Q) the density profile of the trapped Fermi gas is n(r) = n (l — (r /Rtf) 2 ) 3 / 2 , 
and the Thomas- Fermi radius Rtf = Vp/Cl turns out to be the natural length scale in the system. One can also 
introduce the local Fermi momentum Pf(t) = pf(1 — (r /Rtf) 2 ) 1 ' 2 and the density of states on the local Fermi 
surface, vq(t) — mpp(r)/ (2ir 2 ). It should be noted that all these quantities are only slightly influenced by the 
superfluid pairing, because the latter involves only a small fraction of particles (~ T c jep <C 1), with energies close to 
the Fermi energy. 

For describing the phase transition one has to introduce the order parameter which in the case of " singlet" pairing 
is a complex function A(r) ~ (^(r)^/ (r)) e aa >, with e aa > being the antisymmetric tensor. For the "triplet" pairing 
the order parameter is a 3 x 3 complex matrix Ay(r) ~ {o'2^i) a p (i^a( r )9jip '^(r)) , where cr, are the Pauli matrices. 
The time-independent Ginzburg-Landau (GL) equation describes the equilibrium behavior of the order parameter 
A(r) below T c , assuming T c — T <C T c . The critical temperature T c can be found as temperature below which this 
equation has a nontrivial solution. We present the derivation of the GL equation for the trapped Fermi gas in the 
case of " singlet" pairing (g < 0) , relying on the assumption that 

T c ,rj°> > n. 

The derivation of the GL equation for the "triplet" pairing can be performed along the same lines, and should be 
based on the results of Ref. | H]] . 

For the s-wave "singlet" pairing the equilibrium GL free energy (A-dependent part of the free energy) can be 
written in the form 




[A] = Jj~nr Tln/lVexp j - J dr J^(iP + (t,t)i(}-(t,t)A*(t) + h.c.) 1\ , (2) 



F G l[ 

where we use the Matsubara representation, and the symbol < ... >o stands for the average over the states of the 
free-particle Hamiltonian (first term in Eq. (|l|)). In the vicinity of the phase transition the quantity A is small, 
and the second term in Eq.(^|) can be expanded in powers of A. As usual, we perform this expansion up to the 
fourth power. The coefficients of the expansions (kernels) can be expressed in terms of Green function of the normal 
state (without pairing), gL°' (ri, r^). As will be justified below, the order parameter varies on a distance scale Za 
which is much larger than the characteristic distance scale ~ vf/T c of these kernels. Accordingly, Fgl[A] can be 
represented as (hereinafter we use Rtf as a unit of length) 



F GL [A] = R^p [ V-Pr- K o } ( R )I A I 2 " K[ 2) (R) (Ad 2 A* + A*d 2 A - 28 t A* d z A) + K^\A\ 
Jr { \9\ 



(3) 



where R = (ri +r2)/2, r — ri — r2, A = A(R), and 

16tt 2 T 2 ' 



tf«(R) =„ (R)-I%, (4) 



K^ 2 \R) = R 3 TF TJ2 I G ( °i(Rr)G(°)(R,r), (5) 
K[ 2) (R)S l3 = i^irfE J/i r J G -l (R,r)GW(R,r). (6) 
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Here £(x) is the Riemann zeta-function, and the summation is performed over the Matsubara frequencies uj = 
irT(2n + 1), n = 0, ±1, .... The condition Q allows us to use the quasiclassical expression for the product of 

two Green functions: 



m 



\2-KrRTF J 



exp 



2%/2 



V /(l-i?2) 2 + (c/ £F ) 2 + 



l-R 2 



1/2 



(7) 



which can be obtained from the corresponding expression for the spatially homogeneous case, with the replacement 
Pf -> Pf(R)- The validity of Eq. (0) requires the condition (1 - R 2 ) > {fl/T) 2 . The use of Eq. (0) for calculating 
the kernels and Kf is justified by the fact that the pairing takes place only in the central region of the gas cloud, 



and the characteristic size of this region Za <C 1. The main contribution to K\ 
and a straightforward calculation yields 



(2) 
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comes from frequencies |cj| -C £f. 



(8) 



where k = ^7^(3) /48n 2 (VL/T) = 0.13 • (fi/T) < 1. 

The calculation of JTq 2 ' ) is more subtle, because the frequency sum in Eq. (^J) diverges. The divergency can be 
eliminated in a standard way by renormalization of the bare interaction g, and finally one has 



X-K. 



«(R)=i T - w (fl)l„5fSffl 



131 " ' • \a\ T 
Then, the final expression for the GL free energy can be written in the form 



[A] 



R? 



TF 



•)+ 



MR) 



i 



In 



Tc (0) e F {R) 
T e F 



7C(3) 1 

167T 2 T 2 



(9) 



where vq(R) = Uoyl — R 2 , Ef(R) — £f ■ (1 — R 2 ), and the symbol (...) stands for the same combination of 
derivatives of A as in Eq. (||). The analogous expression for the "triplet" pairing can be obtained from Eq. (Q), with 
the replacements A -> T w A 2 /13, |A| 2 -> tr (A+A)^and |A| 4 -> tr (A+A) 2 . 

Since only small distances R ~ y^t <C 1 in Eq. (M) are important for pairing, we will make an expansion in powers 
of R and retain only the largest (quadratic) terms. Then the minimization of Eq. (0) with respect to A* gives the 
GL equation: 



V<9 



2A 



R 1 - In 
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7C(3) |A| : 
8tt 2 T 2 



■A = 0. 



(10) 



We stress once more that Eq.([lO|) is valid under the condition A/T <C 1, which, in turn, implies that T c — T <C T c . 

It is interesting to emphasize that Eq. ( |Io| ) for A is formally equivalent to the non-linear Schrodinger equation for 
the condensate wave function "to in a Bose gas of neutral particles of "mass" 1/2k 2 in a harmonic confining potential 
with "frequency" ("level spacing") 2k = 2k(1 + 1/2A) 1 / 2 . The last (non-linear) term in the l.h.s. plays a role of 

repulsive interparticle interaction, and the third term ln(Tc°^ /T) the role of the chemical potential. Accordingly, the 
calculation of the shape of A is similar to the calculation of the shape of "Jo in a trapped Bose gas, performed in e.g. 
PJQ. But there is an important difference. In the Bose gas the amplitude of is determined by the normalization 
condition. Together with the Schrodinger equation this condition gives the chemical potential as a function of the 
particle number. Hence, for small interparticle interaction the non-linear term is not important at all. In the Fermi 
gas the amplitude of A is always determined by the non-linear term. 

The critical temperature T c for the trapped Fermi gas is the maximum temperature T at which Eq.(]l^) has a 
nontrivial solution. As usual, A — > for T — > T c , and for finding T c the non-linear term in Eq. (|l^) can be omitted. 
Then the GL equation becomes similar to the Schrodinger equation for spherically symmetrical oscillator, and we 
obtain 



T, 



(0) 



T, 



(o) 



In- 



T, 



(0) 



3k« 1. 



(11) 
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One can see from Eq. (|11|) that the critical temperature T c for the trapped gas is only slightly lower than T< 
homogeneous gas with density n . 



(0) 



for the 



As well as in the case of a trapped Bose condensate, the shape of the order parameter in Eq.(lO) is predetermined 
by the ratio z of the non-linear term |A| 2 /T 2 to the "level spacing" 2/5. Since the non- linear term is of order the 

difference between the "chemical potential" ln(Tj ■ /T) and its minimum value (pd|), we have 



2k It 



1/2A) 



-1/2 



(12) 



where ST = T c — T. 

For zCl and, hence, T very close to T c , the non- linear term in Eq.([l0|) does not influence the shape of the order 
parameter, and the latter takes the form of a Gaussian: 



A C (R) ~ cp (R) = (nil) 



-3/4 



exp (-R 2 /2l 



A/ 



(13) 



The linear size of the spatial region where the pairing takes place, Ia — k/vk <C 1, is finite for T — > T c . Moreover, 
Ia ^> £k ~ Q/T c , which justifies the gradient expansion in Eq. (|3|). For finding the amplitude of the order parameter, 
A(i? = 0), in the limiting case z< lwe write A in the form A (R,T) = a (T) (ipo (R) + S(p(R, T)), where ip = ipo + Sep 



obeys the normalization condition J R \ip (R, T)\ =1, and Sip — > for T —>T C . Then Eq. (10) is transformed to 

T + 2A N 



2A 



^ + in(r c /T)(^ + M = « 27C(3)( ^ + ^ )3 



8tt 2 



T, 2 



(14) 



Eq. (|14| ) gives a and 5^ as series of rational powers of ln(T c /T) m ST/T C , the small parameter of expansion being z. 
Multiplying both sides of this equation by ipo (R), integrating over R, and omitting the terms containing Stp, to the 
leading order we obtain 



1 c l 'A 



167r 3 % /2^ 1 T c 
7C(3) 11 T 7 



1/2 



and, hence, 



A(R,T)«A (R,T)=T C 



167r 2 V2 1 T c 
7C(3) n I 7 



^1 
2l\ 



5.15T C 



T r -T 



T c 



-exp 
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(15) 



As well as in the spatially homogeneous case, we have A ~ \JT C — T for T — > T c . For the "triplet" pairing one will 
have Ay(R, T) = Sjj Ap (R,T). It is important to mention, however, that the corrections to Eq.(Q5h, which can be 



obtained from Eq. (|14j), are physically meaningless. They have the same order of mag nitude (oc [8T/T c f/ 2 ) as the 
corrections originating from terms with higher powers of A (for example, A 5 /T*) or higher derivatives, neglected in 
deriving Eq. (|10|). For the same reason one should not go beyond the first term in expanding ln(T c /T) in powers of 
ST/T C . In Fig.l we present the approximate solution Ao (solid lines) and the corresponding numerical solutions of 
Eq. @ (dashed lines) for T c (0) /^ = 5, A = 0.3, and ST/T C = 0.001, 0.01, 0.03. For these values of T c (0) /^ and A 
we have k — 4.4 • 10~ 2 , and Eq. (11) gives the critical temperature T c — 0.87Tc QS> which is only 1% higher than T c 
following from the exact numerica 



solution of Eq.(10) 
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FIG. 1. The order parameter versus R for various temperatures. The solid lines correspond to Aq(-R, T) (15), and the dashed 
lines to numerical solutions of Eq. (10). 

For lower temperature, where z> 1 (but still much smaller than as required by the condition ST <C T c ), one 
can neglect the Laplacian term in Eq. (|10|), and write the approximate solution for the order parameter in the form: 

A (R,T) = T c y^ln^(l - R'/R^f 2 « ZMT C ^EI (l - R^/R^ 2 (16) 

for R<R C = ^J(5T/T C ){1 + 1/2X)' 1 = Z A V2i < 1, and zero otherwise. The solution (|l^) is completely analogous 
to that for the Bosc condensate wave function in the quasiclassical (Thomas- Fermi) regime |l3|,|l4| . 

Eqs. ( |l5|) and ( |l6| ) show that in the vicinity of the phase transition the superfluid pairing takes place only in a small 
central region of the gas sample. This, together with the fact that the superfluid pairing involves only a small fraction 
(~ T c /ep <C 1) of atoms, makes it very difficult to detect the presence of pairing through the measurement of the 
gas density profile. On the other hand, as well as in the spatially homogeneous case, the pairing should influence the 
spectrum of elementary excitations. In this respect we believe that the measurement of eigenfrequencies of oscillations 
of the gas cloud can be one of the most promising ways of identifying the phase transition in trapped Fermi gases. 

We acknowledge very fruitful discussion with G.V. Shlyapnikov. This work was supported by the Russian Founda- 
tion for Basic Studies (grant 97-02-16532). 
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